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Abstract 

Using the method of on-shell recursion relations we compute tree level amplitudes 
including D-dimensional scalars and fermions. These tree level amplitudes are needed 
for calculations of one-loop amplitudes in QCD involving external quarks and gluons. 



1 Introduction and Summary 



The last two years have seen the development of new and surprising techniques for performing 
perturbative calculations in gauge theories. Following the seminal paper by Witten pQ 
(drawing on earlier insights by Nair [2J, quantifying and generalizing the simplicity of some 
tree level amplitudes jSj, the initial efforts have been focused on the relation to twistor 
space (for a review see [4 ). These investigations resulted in a beautiful effective Feynman 
diagram technique [Sj, the so-called MHV diagrams or CSW rules, which can be applied to 
calculations of tree level [Hj, one-loop [Zj amplitudes and amplitudes with external massive 
sources |S]. Recently it has been extended to gravity amplitudes as well [9J. 

Though the CSW diagrammatic method proved much more efficient in calculating am- 
plitudes than the traditional method of Feynman diagrams, they still have the same flavor, 
and it turns out there exists an even more efficient method, that of the on-shell recursion 
relation. This method was suggested in [TU] (see also [12]), and was proven by [UJ. It seems 
applicable to calculations of tree level amplitude (or more generally rational functions) in a 
wide range of theories, and has the flavor of the analytic S-matrix theory, in that it does not 
make use of off-shell structures. As such it seems like a genuinely new way of performing 
calculations in perturbative quantum field theories. 

The recursion relations are reviewed in section 2 below, together with other background 
material. They have been utilized to calculate tree level amplitudes in gauge theories [T3j . 
gravity , amplitudes including massive sources [To"! ITol IT7] and rational functions appear- 
ing in one- loop amplitudes fH] • In all these cases the results obtained are either new, or are 
a more compact form of previously calculated results. The relation to the CSW method was 
discussed recently in [TU] . 

A more challenging task, and one of relevance to upcoming experiments at the LHC, is 
the calculation of one-loop amplitudes. The main ingredient used in calculating one loop 
amplitude is that of unitarity: the multi-valued part of the loop amplitude is determined 
by the tree level results (see for example [2H1E])- The compactness of the results for the 
tree level amplitudes is extremely useful when they are used as input for the calculation of 
one-loop amplitudes. 

The method of generalized unitarity is one of the most efficient general methods of using 
the knowledge of tree level amplitudes in calculating the one-loop amplitudes. It has been 
discussed recently in [221 E3] • In particular then discussion in [22] has concentrated on cut- 
constructible amplitudes (in the sense of [21] )• The more general amplitude has rational 
pieces, which can sometimes be determined separately using recursion relations |25j . 

However, as explained in [2S], a systematic method to obtain the complete amplitude, 
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including any rational parts, is using generalized unitarity in D- dimensions. In continuing 
away from four dimensions, in dimensional regularization, the rational pieces acquire cuts 
as well, and therefore can be constructed using generalized unitarity. The form of the one- 
loop amplitude thus constructed is expected to be compact, reflecting the simplicity of their 
building blocks, the tree level diagrams. 

Motivated by this line of development, we calculate below all the tree level amplitudes 
needed for calculating one-loop amplitudes with up to five partons. These tree level am- 
plitudes differ from the ones previously calculated in that some of the external legs are 
continued to D-dimensions. The one loop amplitudes with 5 partons were previously cal- 
culated in [2E], and we can check the generalized unitarity method against those explicit 
results. Calculations of one loop amplitudes in QCD, with up to six partons, are in progress 

123 ■ 

The outline of this paper is as follows: in section 2 we motivate the set of tree level 
amplitudes constructed here, as the building blocks for the aforementioned one-loop ampli- 
tudes. The method we use is the extension of the BCF recursion relations, and we explain 
the new issues arising when including D-dimensional scalar and fermions. In section 4 we 
exemplify the method by calculating the four point amplitudes in detail. Section 5 is devoted 
to calculations of the five point amplitudes and some checks on them. 

2 Preliminaries 

2.1 Notations 

A massless momentum in four dimension can be written as a product of two (bosonic) spinors 
which we denote by \ a and A Q , so that p adl = a^aPn = ^a^a, where cr M = (1, a), and a are 
the Pauli matrices. We also denote alternatively X a = |A) and A a = [A|, so that p = |A)[A|. 
In the case of several momenta, we also shorten |Aj) and [A< | to |z) and [z| respectively. 

For four dimensional fermions we denote by u±(k) two of the positive energy solutions 
of the massless Dirac equation ftu±(k) = of helicities ±|. These solutions are denoted 
u + {k) — | A) = A Q and u_(k) = |A] = \ a , and they are eigenspinors of 7 5 with eigenvalues =F 
respectively. 

The internal products of these spinors are defined as 

(ij) = u^{ki)u + {kj) [ij] = u + {ki)u-{kj) (2.1) 

which also defines the dual (bra) of each spinor (ket). With this notations |z] = Af and 
[i\ = (\i)a are positive chirality spinors of opposite (=f) helicities, whereas \i) = (A») a and 
(i\ = A" are negative chirality spinors of opposite (±) helicities (note that the internal 
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products are then Lorentz scalars). These consist of the four positive energy solutions of the 
massless Dirac equation. One has the following identities 

= ui + k/i = uj-k/i (2.2) 

where u>± = ~(1 ± 7 5 ) are projections onto the positive or negative chirality subspaces. 



Wave-functions for external gluons can be written in this basis as bi-spinors e a6l 

where q is an arbitrary reference momentum, changing it amounts to a gauge transformation 
on the gluon polarization vector. Similar expression hold for the conjugate part of the 
polarization (e 1 * 1 )"". Working with e a6l amounts to concentrating on a negative helicity Weyl 
spinor, which does not mix with the positive helicity one in a purely massless theory. 

In performing helicity amplitude calculations, it is customary to include the wave-functions 
for internal fermions or gluons in the interaction vertices rather than in the propagator. 
Therefore the propagator is always the scalar propagator \, and numerators which usually 
accompany fermion propagators come about from the two interaction vertices connected by 
the given propagator 1 . 

For other conventions and notation used in helicity amplitude calculations, we refer the 
reader to the review (2l1j. 

2.2 One-Loop Amplitudes with Quarks and Gluons 

The tree level amplitudes calculated here are to be used as basic building blocks in evaluating 
one-loop diagrams in pure QCD, with quark-antiquark pair, and a number of external gluons. 
We briefly describe here this calculation, using generalized unitarity in D = 4 — 2e dimension, 
leaving the details to future work [27 . This serves as a motivation for the particular set of 
tree level amplitudes we evaluate in this paper. 

We will concentrate on color-ordered (or partial) amplitudes, stripping the color indices 
off the external legs. We will therefore not distinguish between various matter representa- 
tions 2 . For the purposes of QCD, the fermions are in the fundamental representation, and 
discussion of the color ordered amplitudes in this context can be found e.g in 

Another standard tool is the supersymmetric decomposition of amplitudes. The one-loop 

amplitudes are easier to calculate if a complete supersymmetric multiplet runs in the loop, 

1 For fermions in a complex representation of the gauge group there is a distinction between fermion 
and anti-fcrmion. As the formalism treats them uniformly, one has to add by hand a minus sign for each 
anti-fermion line. 

2 For example, we will not assume the fermions to be adjacent, as could be assumed for quarks, to include 
the possibility of external gluinos. 
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as loop amplitudes are cut constructible in those cases [24 . This allows us to trade some 
combinations of particles running in the loop for others. For our case, it is sufficient then to 
calculate loop amplitudes with scalars and fermions (but no gluons) appearing in the loop, 
and (as mentioned above) external quarks and gluons. 

The one loop amplitudes can be calculated if we know their singularities in D = 4 — 2e 
dimensions. This requires knowledge of tree amplitudes with two external legs continues to 
D dimensions. Using four dimensional helicity regularization (i.e, where the momenta but 
not the polarizations are continued to D dimensions), the D dimensional external momenta 
can be thought of as 4 dimensional massive momenta [2H! • Inspection of the quadruple and 
triple cuts of the aforementioned one-loop amplitudes leads to the following 4 sets of tree 
level amplitudes (see figure 1), which will be calculated below: 

• Amplitudes including 2 massive (or massless in D = 4 — 2e dimensions) scalars, and 
some number of gluons. These were already calculated in the first reference in |15j . 

• Amplitudes including 2 massive fermions (which we denote by A), and some external 
gluons. We will label these type-A amplitudes. 

• Amplitudes including 2 massive scalars, and massless gluons, accompanied by an ex- 
ternal quark- ant iquark pair. We will label these type-B amplitudes. 

• Amplitudes including a massive scalar and massive fermions, with external massless 
fermion and a few gluons. Those amplitudes will be labeled type-C. 



Figure 1: The relevant tree diagrams, with varying number of gluons, are ordered from left 
to right. The massive legs are denoted in boldface lines, solid lines denote fermions, dashed 
lines scalars, and wiggly lines gluons. 

We note that all amplitudes involve two massive legs which are adjacent. 

In an effort to use a uniform notation, we will denote the momenta of external 3 fermions 

by k, external scalars by /, and external gluon by p. In addition, as explained in the next 

subsection we need to distinguish between four dimensional momenta (denoted by small 

letters) and D-dimensional ones (denoted by capital letters) . 

3 We will not use uniform notation for intermediate state momenta that are eliminated from the final 
result anyhow. 
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2.3 D-dimensional Fields 



We find that in order to use D-dimensional unitarity, we have to use helicity methods in 
calculating tree level amplitudes in D-dimensions. In the FDH regularization scheme, the 
momentum is continued to D-dimension. Every D-dimensional momentum P can be de- 
composed as P = p + fi, where p is the four dimensional component, and fi is a compo- 
nent in a formal (— 2e)-dimensional orthogonal space 4 . Working in mostly minus signature, 
P 2 = p 2 — fi 2 , so on shell massless momentum (P 2 = 0) is equivalent to four- dimensional 
massive momentum p 2 = /i 2 . Therefore for scalars, working away from 4 dimensions is 
equivalent to adding mass to the scalar field. 

Each loop momentum integration can be decomposed as ?f~p7 — (f^ T^r) 3 ^ ' so ^ e mass 
fi is always integrated over, and the e dependence of the amplitude is generated from the 
\i dependence of the integrand. Similarly, each vertex is accompanied by a delta function 
imposing momentum conservation. As all momenta are now D-dimensional, those are also 
D-dimensional delta functions. It is therefore necessary (in general) to regard the momenta of 
internal lines as being D-dimensional as well; for tree level amplitudes they are just given as 
linear combinations of external momenta, given by imposing all the momentum conservation 
constraints. 

For internal fermionic lines one always has to sum over the intermediate spinor wave- 
functions, so choice of basis is not necessary. We will therefore use the notation \P}, {P\ 
to refer collectively to these wavefunctions. The sum over the intermediate wavefunction is 
performed using the identity 



Note that now If has one component (^) that preserves helicity and one that flips it (/I), 
similarly, in D-dimensions the components p 1 and i/l behave differently with respect to chirality, 
{A 7 5 } — whereas [/Jl, 7 5 ] = 0. 

To mimic the four dimensional helicity methods, we want to utilize helicity-like states 
for external fermions, even when they are D-dimensional. As we keep 7 5 four dimensional, 
we can still use chiral basis which we denote as before by \P) = lj + \P} and \P] = uo-\P} 
(and similarly the conjugates {P\ and [P\). However the states of definite helicity, \P) and 
[P\ (or similarly \P] and (P) now mix with each other. Indeed, the basis vectors \P) and 
\P] do not individually satisfy the massless Dirac equation in D-dimensions. That equation 
written in terms of these Weyl fermions is, 



4 We will use the capital letters for D-dimensional momenta, and small letters for their four dimensional 
components. The (— 2e) component of all D-dimcnsional momenta is always ±/i. 



P}{P\ = t? 



(2.4) 



i>\p) + M = o 



i>\p] + i/i\p) = o 



(2.5) 
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which is consistent with the mass-shell condition p 2 = fi 2 . 

Nevertheless one can assemble the physical amplitudes (with external wave functions 
\P} = \P) + \P]) from the ones calculated here, as helicity violation is limited to insertions 
of i/i in fermion lines. In the course of using the recursion relations we demonstrate this 
process, which is also relevant for the calculation of one-loop amplitudes using generalized 
unitarity [21] • Whenever one encounters an intermediate D-dimensional fermion, one can 
write the numerator of the propagator as 

P- = \P}{P\ = (\P) + \P])((P\ + [P\) (2.6) 

and make use of the partial amplitudes with helicity states, the ones we calculate here. 

Note however that the propagator has both helicity preserving and helicity flipping parts. 
The helicity preserving parts are the usual propagators, usually drawn as connecting ± states 
to =F states, 

\P)[P\ = u+j> \P\{P\ = u>-j (2.7) 

whereas the helicity flipping parts are new, and are the part of the propagator that connects 
± states to ± at the other end of the propagator. They arise from the identities 

\P)(P\ = u-t \P][P\=lu + ^ (2.8) 

The main advantage of using chiral external states is that one can use the simple expres- 
sions for the gluon polarizations, 

m = ^^(\Q)[k\ + \k](<l\) = ^ + )a a + (e + r 

frw = -^^ ] (\m + m\) = (t-u + (e + r (2.9) 

This polarization contracts into the fermionic states that accompany the gluon in an inter- 
action vertex. In case at least one of these states is a Weyl fermion, one of the terms in the 
polarization vanishes (which one depends on the chirality of the fermionic state). Note that 
in this case the chirality of the other fermion is determined, even if it is a D-dimensional 
(and thus Dirac) fermion. 

To summarize the step of the calculations, we will be using the BCF recursion relation to 
find an expression for helicity amplitude which is valid in general dimension D (in the FDH 
scheme). Subsequent manipulations will depend on which external (and internal) legs are 
taken to be D-dimensional, those include translating the helicity states into momenta, using 
mass shell conditions and trace identities to simplify the results. We will be very explicit 
in calculations of the four point amplitudes in section 3, to demonstrate the issues involved, 
and will be less detailed in deriving the five point amplitudes. 
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2.4 Recursion Relations 



To evaluate the tree level amplitudes, we will use the recursion relations first discussed in 
[Til] and proven in ^JJ. The proof utilizes the analytic properties of rational functions, and 
can be then generalized to massive particles [2H|, and to purely rational loop amplitudes [7j, 
or to calcualting the rational parts of loop amplitudes JH] • Here we slightly generalize it for 
the case of D-dimensional fields. We briefly review the method as needed for our purposes, 
highlighting the slight differences arising in our case. We refer the reader to a more detailed 
discussion in the original papers. 

We will be discussing tree level amplitudes A n (p\, ■■■,p n ) of n external on-shell particles, 
n — 2 of which are massless in 4 dimensions, or two are taken to be massless in D = 4 — 2e 
dimensions (or equivalently massive in 4 dimensions). The recursion relations depend on 
choosing two of the external momenta (labeled i,j) and "marking" them. 

Now, define a function A(z) to be the amplitude evaluated at the shifted momenta 5 
Pi = Pi + zrj and pj = pj — zr], where r\ is a null vector orthogonal to both pi, pj. This ensures 
the same mass-shell condition applies to the shifted momenta. The shifted momenta are now 
null and complex. We will always choose at least one of the marked momenta to be massless 
in 4 dimensions, as this simplifies the analysis, and is sufficient for our purposes 6 . In case 
both marked momenta are null, we can write them as product of spinors: Pi = AjAj and 
Pj = XjXj, then rj = XjXi. In this case the shift amount to shifting the spinors Aj — ► Aj + z\j 
and Aj — > Aj — z\i, leaving Aj and Aj intact. In case that pi is massive and pj = \j)[j\ is 
massless, we have rj = \j\pi\j) (and similarly for i j). 

One then divides the n external momenta to two cyclically ordered groups, which are 
labeled L = {p r , ...p iy ...,p s }, R = {p s +i, ■■■Pj, ••-,Pr-i}- As is indicated the groupings is such 
that pi e L and pj e R, and we will sum over all such groupings. We denote by p = p r + ... +p s , 
the momentum flowing in the channel between the L, R groups of momenta 7 , and p = p+zrj is 
the shifted intermediate momentum. The shift variable z is chosen to impose the appropriate 
mass-shell condition on the shifted intermediate momentum p. For uniformity of notation, 
we impose the same mass-shell condition imposed on external legs: p 2 = for purely four 
dimensional momenta, and p 2 = fi 2 for components of D-dimensional momentum P. 

5 Note that the momentum conservation constraint is unaltered by the shift, so A(z) can be calculated 
using perturbation theory. 

6 For simplicity we take rj to be purely four dimensional null momentum, even when D-dimensional mo- 
menta are involved. 

To conform with the notation in this paper, we use p for four dimensional intermediate momentum, and 
P for a D-dimensional such momentum. 
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The BCF recursion relation is then 

An(pt, ...p n ) = ^2J2 A l(Pt, -Pi, -,Ps, -p h )\A R (p~ h ,p s+1 , ...,pj, ...,p r ~i) (2.10) 

L,R h P 

The first sum is over all possible groupings of external momenta, as described above. The 
second sum is over all possible intermediate states (depending on the matter content of the 
theory), and their helicities 8 h. The amplitudes Al, Ar depend on shifted momenta as indi- 
cated, and are momentum conserving on-shell amplitudes, albeit with complex momentum; 
they include an additional external leg with momentum ±p and the appropriate helicity ±h. 
Note also that the momentum p appearing in the propagator is unshifted. 

The validity of the recursion relations depends crucially on one technical assumption, 
that of the vanishing of the function A(z) as z — ■> oo. This depends on the helicities of the 
marked momenta In general, the helicities (hi,hj) must chosen so that hi > hj, with 
additional constraints when choosing quarks. In particular, for two gluons we cannot choose 
(— , +) [ID]; when choosing a gluon and a scalar, positive helicity gluons must be particle i, 
while negative helicity gluons must be in the position j ^3] . When "marking" a quark, we 
cannot also choose an adjacent quark nor an adjacent scalar, and for an adjacent gluon-quark 
pair, we have similar rule as above: positive and negative helicity gluons must be chosen as 
i and j, respectively [Tfi] . 

We would now like to resolve the issue of choosing an adjacent gluon and quark (in that 
order) with helicities (+, +). While claims that this choice is invalid, we will argue that 
it is in fact allowed, as stated in We will consider A(z) as a sum of Feynman diagrams 
where the momenta of the gluon % and adjacent quark j depend on z, and we will write 
pi = p(z) and pj = k(z). When i and j share a vertex, the z-dependance of those diagrams 
are completely determined by the corresponding polarization vector and wavefunction. This 
class of diagrams contributes the factor [k(z) \ (e + (p(z))). We recall that [k(z)\ oc z and 
e + (p(z)) oc l/z, and so naively this term would not decay for large z. However, we point 
out that 



so this type of diagram does vanish for large values of z. When the particles % and j are 

separated by any single (^-dependent) propagator, it is straightforward to check that similar 

cancelations occur in the numerator whenever an amplitude poses the threat of not decaying 

at large z. Finally, any further lengthening of the "z-path" is harmless since the only 

components which grow with z are the scalar-gluon and the triple-gluon vertices. However, 

8 For D-dimensional momenta helicity is not well-defined, but one still have to sum over an appropriately 
chosen basis, as described above. 




(2.11) 
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introducing such interactions must be accompanied by the appropriate bosonic propagator, 
which then compensates for the z behaviour of the vertex. Thus, we find that the choice 
(+, +) is valid for an adjacent gluon-fermion pair (in that order), and similarly for the choice 



The new element in our calculation is the continuation away from four dimensions. How- 
ever, the continuation to D dimensions does not affect the analytic properties of A(z) or 
any of its ingredients, so the the recursion relations remain valid, as do the choices for 
given in the literature. This fact was already noted in |16] . We point out that the only 
qualitatively new ingredient in D dimensions is the helicity-flipping fermionic propagator. 
However, this in fact has even better large z behaviour then its standard helicity-preserving 
counterpart, since the numerator fi is constant, as opposed to If which is linear in z. 

To summarize, in the following we will use the BCF recursion relations, always marking 
momenta in configurations that are proven to be allowed (that is, when the vanishing of the 
boundary term has been established, as described above). 

2.5 Primitive Vertices 

We are interested in calculating tree level diagrams involving massless quarks and gluons, 
and massive scalars and fermions. Using the BCF recursion relations ^0], we need as basic 
building blocks a few cubic amplitudes. We list here the vertices we need for our calculation. 
In every case these are simply constructed from the gauge theory action by contraction with 
the external wavefunction of the appropriate helicity. In all the cases listed below we treat all 
momenta as incoming. All vertices are written for k) cyclically ordered, the expressions 
for the other cyclic ordering k) differ by an overall sign if are fermions. 
The primitive vertices we need are: 

• gluon-fermion vertex (figure 2): 



where q ia an arbitrary reference null vector with q = \q)[q\- The vertex is independent 
of this choice as a consequence of gauge invariance. As each one of the primitive vertices 
is gauge invariant, these reference vectors can be chosen independently for each vertex. 

• scalar-fermion vertex (figure 3): 

negative helicity external legs: (ij) 




negative helicity gluon k: 
positive helicity gluon k: 



{ik)[qj} 
[qk] 



jjg) [Ml 

(qk) 



(2.12) 
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Figure 2: Gluon-Fermion primitive vertex, the fermion helicities are drawn, the gluon can 
have positive or negative helicity. 



positive helicity external legs: —[ij] 



(2.13) 




Figure 3: Scalar-Fermion primitive vertex, The 3 external legs have the same helicities. 



scalar-gluon vertex (figure 4): 

negative helicity gluon k: 
positive helicity gluon k: 



{JAM 

[qk] 

(aim 

{qk) 



(2.14) 



gluon cubic vertex (figure 5): 
MHV vertex: ^ 

MW vertex: ^ (2.15) 
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Figure 4: Gluon-scalar primitive vertex, the scalar helicities are drawn, the gluon can have 
positive or negative helicity. 





Figure 5: Gluon MHV (left) and MHV (right) cubic vertices. 



There are other non-vanishing cubic vertices which we will not need, therefore we will not 
present them here. Additionally, for D-dimensional fermions one utilizes the full expression 
for the gluon polarization, resulting in additional terms in the interaction vertices which we 
write down when we use them below. 

2.6 Checks on the Amplitudes 

Though the on-shell recursion relations are proven, it is still useful to perform a few checks 
on the resulting expression, verifying the various intermediate steps leading to the final 
expression. These steps include choosing marked momenta (such that A(z) vanishes as 
z — > oo), choices of various reference momenta, and straightforward (but sometimes tedious) 
algebra. 

The first check one can perform is comparison with the result of Feynman diagram 
calculation. We have checked our expressions against such calculations for all the four point 
amplitudes and some of the five point ones. Typically the recursion relations yield much 
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more compact expressions for the amplitudes, so the main complication is to reduce the 
complex Feynman diagram result to the simpler expression. 

Another check we have performed is (some of) the collinear limits of the amplitudes. The 
collinear limits are a subset of the multi-particle poles which occur at tree level amplitudes. 
As the sum of two neighboring momenta becomes on-shell the amplitude factorizes in the 
appropriate channel. All our amplitudes have poles at the right location, and in some cases 
we have checked explicitly that the residue of the pole is the expected one. We exemplify 
the collinear limit in the appendix. 

Finally, in the limit fj, — > all external legs are four dimensional. In some cases the 
amplitudes are known in that limit, and we reproduce these results. 



3 Four Point Amplitudes 

The amplitudes with four external legs are fairly simple, and can be checked explicitly against 
Feynman diagram calculations. We present the details and the results in this section as a 
demonstration of the technique and the new issues arising when D-dimensional fermions and 
scalars are included. 

In addition, all such amplitudes can be seen to have the correct factorization limits. 
Indeed, to see the singularity structure it is sufficient to inspect the denominators, we easily 
see that they vanish if and only if the sum of two adjacent momenta becomes on-shell 
(so spurious singularities are absent). Checking factorization amount to verifying that the 
residue of these poles is the expected one 9 . 

3.1 Type- A Amplitudes 

These amplitudes have two adjacent massive fermions (of momenta and helicities K^K^) 
and two adjacent gluons (of momenta pi,P2)- We discuss all helicity configurations in turn. 

The first case of where the gluons are of opposite helicities is the amplitude Ai{K^ , K 2 , pf , P2 )• 
We choose the marked momenta to be (i, j) = (pt,P2 ), as this is one of the configurations 
for which there is a general proof of vanishing at infinity. In this case there is only one pos- 
sible diagram appearing in the recursion relation, with intermediate fermion of momentum 
P = K 2 + p\. One gets 

9 The term collinear limit is inaccurate when massive momenta are involved, what we really mean is 
two-particle factorization limits. 
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where q±, q2 are two reference momenta. Note that the helicity of the intermediate states of 
momentum P is determined by that of the external momenta. 

We choose q\ = 2 and q 2 = 1, and use |1] = |1] and |2) = |2), then this becomes 

(21)[12]P 2 (21)[12]P 2 1 ' ' 

This result is valid in D-dimensions, with states such as \K 2 ) are defined in section 2. Note 
also that in this case the relevant part of P appearing in the intermediate state is p, as 
[l|/i|2) = by chirality selection rules. 

The next step we simplify the result by taking the gluons to be massless in four dimen- 
sions, and the fermion momenta to be massless in D-dimensions (with — 2e component p). 
Then the four dimensional component of P is p = p\ + k 2 , and the — 2e component is /i. This 
result can be shown to be identical to the one obtained from using Feynman graphs 

MKtM-) ^^ ,3,8) 

A slight variation of the same calculation yields 

A,(K 1: K 2: 1 ,2 )- p2(12)[21] " { Pl + k 2 y-^ (3 ' 19) 

where the intermediate momentum P = p\ + K 2 in the only diagram contributing to the 
recursion relations. 

Additionally, for these gluon helicities there could be one helicity flipping amplitudes, 
using the same choices of marked and reference momenta one gets: 

^/ f ?4^r,^I<« ,3,0) 

However, using |P)(P| = uj_i/l and (2|^|2) = this amplitude vanishes. 

Now, if the gluons are of the same helicity, we get the amplitudes A^Ki, K 2 ,pf, p 2 ), 
where the massive momenta K\,K 2 are of either chirality. In this case both external states 
are of the form \Ki} = \Ki) + \Kj\. Choosing the marked momenta to be = {pf-,P 2 ) 
gives one possible intermediate (D-dimensional fermionic) state with P = K 2 + p\. 

When contracting with the gluon polarizations / + one gets the interaction vertex 

-\-((Pq)[iK 2 } + [Pi}(qK 2 )) (3.21) 
with a similar expression for the other interaction vertex in the diagram. We get 
A 4 (K U K 2 , l+,2+) = -L- ({P qi )[iK 2 } + [Pi]( qi K 2 )) 4 



{qil) V — - l — "J P 2 

((-Pq 2 ) [2K,] + [-P2] (q.K,)) (3.22) 
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There are four terms, corresponding to the four possible helicity assignments for K±, K 2 . 
As the amplitude vanishes in four dimensions, we know that only helicity flipping terms has 
to be retained. Starting with Ki, K 2 we get (choosing q± = 2 and q 2 = 1), 

A l{ K,,K^2 +) = JlMMmEA (3,3) 
using |P][P| = uj + i/l , and the more standard projections gives 

(21><12> 1 1 

As i/l anti-commutes with / , )2 , those combine to give (1 + 2) 2 — ffi . The last term gives 
vanishing contribution (using e.g. [P|^|ifi) = by momentum conservation), therefore 

A (K- K- 1+ 2+) - [12] {KMKl) (3 25) 

' K * ' 1 ' 2 } " J12) (k 2 + Pl y - ^ (3 - 25) 



This matches the result quoted in jzy]. Similarly for the last helicity configuration one 
obtains 

im 2 (3 , 6 ) 

(12) (k 2 + P1) 2 ~ P 
3.2 Type-B Amplitudes 

These amplitudes have two massive scalars (of opposite helicities) with momenta L 1; L 2 , and 
two massless fermions of opposite helicities, and momenta k\ = X\Xi and k 2 = A 2 A 2 . As 
we are interested always in adjacent massive legs, the only non- vanishing helicity preserving 
configurations have cyclic ordering of momenta (kf, k 2 , L±, Lf) or {kf , k 2 , Lf , L 2 ). The 
helicity violating configuration is (kf, kf, Lf, Lf). 

For all these configurations we choose the two marked momenta to be (i,j) = (ki,Li). 
In this case there are two possible grouping of momenta, or two possible diagrams in the 
recursion relation. In one of them the intermediate momentum is P = k\ + L 2 , and the inter- 
mediate state is a D-dimensional fermion. In the other diagram the intermediate momentum 
is q = ki + k 2 , and the intermediate state is a gluon which can be of positive or negative 
helicity. The first set of diagrams can lead to helicity flipping via the D-dimensional inter- 
nal fermion, whereas the diagrams with an internal gluon only lead to helicity conserving 
amplitudes. 

The first amplitude A^(kf , k 2 , L\ , Lf) can be written as a sum of three terms 

[k 1 -P}(Pk 2 ) (k 2 qi)[k\-q] 1 (q\l 2 \q 2 \ 



Ai{k 1 , k 2 , L x , L 2 



P 2 (9i -q) Q 2 M 

[q3q] Q 2 {m) 
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The last term vanishes if we choose g 3 = k\ (note that (kiq) ^ 0). For the middle term we 
choose qi — q 2 — ki, the amplitude is then 

72 > { lh\h\k 2 ) + | (3.27) 

using elementary consideration the numerator can be simplified 

[kiWMlh) = ~2(h ■ k 2 )[k 1 \l&\h 2 ) = -2(h ■ q)[kM\k2) (3.28) 

therefore 

MK,K,L T ,m = Nwj ,, + l! '.-/ jkTTw) < 3 ' 29 > 

which can be easily checked to be the result obtained from Feynman diagrams. 

Using similar reasoning, the the other helicity conserving amplitude evaluates to be 

Mkt,K,Lt,I^) = ^f^ (3-30) 

Note the absence of singularity as k\ + L 2 becomes on -shell, as there is no appropriate 
helicity assignment for the would-be on-shell intermediate state. This result can be easily 
checked to be the result of a sum of two Feynman diagrams. 

As before, the intermediate state is effectively four dimensional for helicity preserving 
external states. We now evaluate the helicity flipping amplitude, for them the only possi- 
ble intermediate state is fermionic and its propagator is effectively ^, which simplifies the 
calculation. The primitive vertex coupling scalars to fermions is unmodified, giving: 

MktktLtLt) jk^^,_^^ (3.3!) 

However, for the purpose those type B amplitudes where the "helicity" of the scalars is 
flipped are not used as ingredient in the one-loop calculation we are ultimately interested in. 
We therefore only consider the case where the two massive scalars have opposite "helicity" . 
This is useful because often this eliminates the fermion helicity-flipping terms as well. 

3.3 Type-C Amplitudes 

There are a few helicity configurations relevant here. We will sketch the calculation and give 
the results. In all cases we take the fermionic momenta to be fci, K 2) the scalar momenta L 
and the gluon momentum p. 
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For A^kf, K 2 , L + ,p + ) we choose the marked momenta to be = (p + , kf). There is 
one possible grouping of momenta for which the intermediate state is a D-dimensional scalar 
of momentum Q = p + L. The recursion relation reads 

A 4 (kt,K+,L + ,p + ) = -<jM±-[k lK2 ] (3.32) 



choosing q — k\ gives 



(3.33) 



(hp) Q 2 

where we already projected onto the appropriate component of $ , namely /. 

The last step involved calculating [X2&1] which will be used repeatedly below. Using the 
result 

one finally gets 

MMM- J^l^ (3,5, 

For the amplitude A^kf , K 2 , L + , p~) we choose the marked momenta = (kf,p~), 
resulting in similar calculation (where the helicity of the external gluon is flipped). The 
result is 

Mkt , K} , L ^ ^mA_ ] ( ,sa) 

Note that for these two amplitudes, where the fermions are adjacent, the intermediate 
state is massive scalar and consequently there is no helicity-flipping intermediate state. 

For the amplitude A^kf ,p + , K 2 ', L + ) we choose = (p + , kf), therefore the interme- 
diate state is of momentum Q = p + K 2 , and 

A 4 (kf,p\K+,L+) = - [m ^ p ~ ^ ^ IQh] (3.37) 

choosing q = k\ gives after some elementary algebra 

A 4 {ki,p + ,Ki,L + ) = f (3.38) 

(k 1 p)(pK 2 ) 

As the intermediate state in this amplitude is D-dimensional fermion, there is a similar 
amplitude with opposite helicity for one of the fermions, namely 
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after calculating we end up with the amplitude 

M^M^ J™f\^ ,3,0, 

Finally for the amplitude A±{k1[ ,p~, K 2 , L + ) we choose = (ky,p~), giving 

Mkt,P~,K},L + ) = -[k - P]^ {m { ^ 2] ( 3 - 41 ) 
with P = K 2 + p, choosing q — k-y gives 

AAk + v- K + L+) - [hK2] Wlkl] (3 42) 

A^k^p ,K 2 ,L )- ^ ^ + 1)2 _ ^ [6A2) 

Once again, since the intermediate state has helicity flipping part, we could get the 
additional amplitude 

Mkt,P~,K^L + ) = i^- p ]^ {K f^ P] ( 3 - 43 ) 
However, choosing again q = k±, and using [A^l/ll/cx] = 0, this amplitude vanishes. 

4 Five Point Amplitudes 

We list below the results of the calculation of the relevant five point amplitudes, and the 
checks they satisfy. In all cases the calculation follows the lines of the corresponding four 
point amplitude calculation, and we omit the details for brevity. 

4.1 Type- A Amplitudes 

These amplitudes include two adjacent D-dimensional fermions of momenta K\, K 2 ,an.(\. three 
adjacent gluons of momenta Pi,P2,P3- The results are 

a (k- k- v+ v+ v+) - ^MK^mm+Hm 

' 2 ' Pl ' P2 ' P3j " (X 2 1)(12)(23)(3K 1 )[X 1 3] 

There is no singularity as K\ + K 2 becomes null, since no helicity assignment exists for an 
intermediate state. The result is symmetric when exchanging (1,3) and (Ky,K 2 ), therefore 
there are only two collinear limits to check. We perform these checks in the appendix as a 
demonstration. 

Note also that the amplitude vanishes in four dimensions (setting fi = 0) as it should, 
and therefore only helicity-flipping parts exist. The other possible helicity assignment for 
fermions gives 



As(Kt,Kt,p+,pt,pt) 



(X 2 1)(12)(23)(3X 1 )[X 1 3] 
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The rest of the type A amplitudes are given as sum over two diagrams, one with internal 
gluon and one with internal D-dimensional fermion 



51 1 ' 2 ' p i' p 2.Ps; - (x 2 l)(l2)(23)(3X 1 ) [Z\{kl + kfi\K 2 ) 



A 5 (K 1 ,K 2 ,pt,p 2 ,pt) 



A 5 {K+,K+,p+,p 2 ,p+) 



[U^SllK^iK.K,) [3\(kl + kfi\K 2 ) 

(k 2 2)(2k 1 ) 

(2i)(K 2 + Pl y(p 2 + p 3 y(K l +p 3 y 

[i^iK^m 

[UnK.K^K^) [3\(kl + ty)\K 2 ) 



A 5 { K x , K 2 , Pl , p 2 , p 3 ) - (12> (23) ^ ! ( ^ + ^ ! i} + 



fi 2 [12](K 2 3)[2K 



MKt,K+,pt, pt, Ps) 



A 5 (Kf,K 2 ,pt,P 2 ,Ps) 



A 5 (K+,K+,p+,p 2 ,p 3 



(12)[23] (K 2 +p l y(K 1 +pzf 

(lK^jK^kj + k^l) [K 2 \(kj + km 
(K 2 1)(12)(23)(K 1 K 2 )[K 2 K 1 ] [K 2 \{ki + fc^)|3) 

[23] [^ 2 2] 2 (K 2 1> [KMlKz] 



(p 3 + ^i) 2 (^2 +Pi) 2 [AT 2 1 (A;/ + A$|3) 
4.2 Type-B Amplitudes 

These amplitudes have two massive scalars (of opposite helicities) with momenta L\, L 2 , and 
two massless fermions of momenta ki,k 2 , and an additional gluon of momentum p. The 
helicity flipping part of the fermion propagator are less relevant in this set of calculations, 
since fermion helicities are correlated with that of the external scalars. The results are 
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a 5 ( p + , Lr, Lt, kt, k 2 ) = {kMp] [kimh) ( {hmp] 



(hp) {p\U\p} \[h\m + v 2 )\p] (hk 2 ) (h + L 2 )\ 

n 2 {ph] 3 



[ki\(Wl + ^)\p] [k 2 h] (L l + L 2 y 



Mv~ L- L+ k+ k~) ~ ( mih] + 1 ^ x 

x ( \kM\v) \kM\k 2 ) - A^mv) \ + /i 2 (ph) 2 (Ph) 



<pM|p] [p2] ; (p|(^ + fi 2 )\k 2 ) (hk 2 ) {U + L 2 )^ 

(k 2 \U\p\ 2 (ki\U\k 2 ] 



•Az,(p , 1^ , L 2 , k l , k 2 



(hk 2 )(k 2P ) (p\U\p}[ki\Wl + ^)\p] 

A* 2 [pk 2 ] 2 [ph] 

[hk 2 ] [A; 1 |(/i// + /i 2 )b](L 1 + L 2 ) 2 

(p\lj\k 2 ] 2 (hmm 

[kik 2 ][k 2 p](p\U\p]{p\m + » 2 )\ki) 

^ 2 (Phf 

(hk 2 ) (pim + ^lkdiU + Lrf 



^5(^1 , P + , L 2 , k 1 , k 2 ) 



(hp)(pk 2 ) \(k 2 m\ki) (L 1 + L 2 ) 2 J 
1 V 2 [p\U\k 2 ) 2 
(k 2 \m\ki) (L 1 + L 2 y(L 1 + k 2 y 

(k 2 \m\ki) 



(hp)(pk 2 ) (Lx + L 2 ) 2 
We also find that all amplitude where fermion helicity is flipped are vanishing. 
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4.3 Type-C Amplitudes 

These consist of massless fermion of momentum ki, massive fermion of momentum K 2 , 
massive scalar of momentum L and two massless gluons of momenta Pi,p 2 . The results are: 

,2 



A 5 (pf,pt,kt,L + ,K 2 



A 5 (pt : p 2 ,kf : L + : K 2 ) 



(K 2 l)(12)(2k 1 ) 
(K 2 l)(12)(2k 1 )(k 1 + Ly 

(2K 2 )[hm) 2 

(K 2 1){12) [hltf + kfilKt) (h + Ly 

Ahlf 

[12] [2**] [h\(Jt + H)\K 2 ) (fcx + L)2 

(K 2 2) 2 [1|^|A; 1 ](2^|A; 1 ] 
(K 2 l)(12)[k 1 \jl\K 2 )[K 2 l}(k 1 + Ly 

(K 2 \([ + kmiMhBh} 2 



A 5 { Pl ,pt,kt,L + ,K+) 



[iir 2 l][12][2A;i][A;i|(/ + H)\K 2 ){K 2 + Lf 

(ih)(i\t+mK 2 ] 2 



,pt,kti L+ i K 2 ) 



,P 2 ,Ki L+ i K 2) 



A 5 (pt,kf,Pt,L + ,K+) 



A 5 (pt,kf,pt,L + ,K 2 ) 



(12)(2fci) (hlijl + ty)\K 2 ] (K 2 + L) 2 

y 2 [K 2 2f 

[12][X 2 1] [K^ + k^ih + L) 2 

[^2|(/ + ^)|A;i>[iir 2 l][12](A;i + L) 2 

[K 2 h} [h^+jh + LYlK,} 
[12] [K 2 l][2h] (h + L) 2 



(K 2 l)(lk 1 )(k 1 2){2m 

(^ 2 fc 1 )[l|^(fcj/ + /x 2 )|2] 
{K 2 l){lk 1 ){k 1 2){2\^\2\[lK 2 \ 
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(lh) (2\(tlty + ^)\l)(K 2 + Ly 



A 5 (pt,kf,p 2 ,L + ,K 2 ) = 



[^2](K 2 1)(2^|2] (2\(jtK£ + H 2 )\l) 

(K 2 \m+m+^)mmhm\h] 



As(p 1 ,kt,p 2 ,L + ,K 2 +) 



(l|(A;i/ + /U 2 )|2)(lir 2 )[A; 2 l][A; 1 2](2^|2] 

mkj + f)\K2} 2 
[/f 2 l](fc 1 2>(2| / i|2][l|(^ + ^)|2] 

+ [lk l ][2\im + ^)\\]{K 2 + Lf 
[K 2 k 1 f{2^\k l \ 



[tf 2 l][lfci][fci2](2W2] 



+ r _| ,, + ^ [2|(/^ + ^ 2 )l^i](^i|^|^2] 



A(^,p^,p 2 f ,^+,^ 2 



A(^,^,p 2 -,l + ,^ 2 + ) 



A 5 {kf,pt,p 2 ,L + ,K 2 ) 



(A; 1 1)(12)(2^|2][(X 2 + ^ 1 ) 2 -^ 2 ] 

" (A; 1 l)(12)(2^|2][2|(/A;i + / i 2 )|X 2 ](ir 2 l) 

(fc 1 2)[A: 2 |^ + fc^)|2)(2|/|Ar 2 ](A: 2 2) 
" (k 1 l)(\2)(2\{m + ^)\k 1 ){K 2 + Lf 

(A: 1 |A:i|ir 2 ][A: 1 l](2^|l] 2 
[12][2^|2)(A; 1 |(^ + / x 2 )|2)(A; 1 + K 2 ) 2 

(if 2 fci)[l|^N(2^|l] 3 

" [X 2 l][12][2^|2)(2^(/ + ^)|^ 2 )(2|(/A;i + / x 2 )|A; 1 ) 

(im 2 [kiK 2 ] 2 (lK 2 ) 

1 ; (12)(2i/|2][A: 1 |(^ + /i 2 )|2][(A; 1 + K 2 ) 2 -/i 2 ] 

mk^^)\K 2 \[2k x f 



+ 
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A 5 (kf,p 1 ,p^,L + ,K 2 



A 5 (kt,pt,pt,L + ,K+ 



[^il][12][2|/|2)[(/ci + K 2 ) 2 - fi 2 } 

To summarize, we have applied the BCF recursion relations to amplitudes which involve 
D-dimensional fermions and scalars. These D-dimensional particles behave in most respects 
as their massive counterparts in four dimensions. We have used this formalism to obtain four 
and five point amplitudes where two of the particles have been continued away from four 
dimensions and the remaining particles are (on-shell) gluons and massless fermions. Our 
results posses expected factorization properties and have passed other consistency checks, as 
outlined in section 2.6. The tree-level amplitudes we have computed here are the building 
blocks (to be assembled using generalized unitarity in D dimensions) of the rational (non- 
super symmetric) terms of the one-loop amplitude of two (massless) fermions and up to three 
gluons. We leave this task for future work [2*7] . 

Appendix 

Let us discuss the factorization limits of the amplitude 

A 5 {K,,K 2 , Pl ,p 2 ,p 3 )- {K2l){u){23){3Ki)[Ki3] 

As mentioned in the text, there is no singularity in the channel where K\ + K 2 becomes 
on shell. In addition there is a symmetry of exchanging (1, 3) and (Ki, K 2 )- This leaves two 
channels to check, when p\ +p 2 becomes null, or when p% + K\ becomes null in D-dimensions 
(so it approaches the mass shell condition (p 3 + k\) 2 = fi 2 in four dimensions). 

For the first limit we denote p — Pi + p 2 , the amplitude as p becomes light-like factorizes 

to 

1 



A 5 (K 1 , K 2 ,pj ,pj ,pj) -»• A 3 (pJ ,P 2 ~P ) 13 Ai(p + ,pt, K 1 , K 2 



p 

The four point amplitude is of type A, so we get: 

1 [lp][p 2 ] (12>[21] (p3}(K 1 +p 3 y 

To get to the right form we multiply both numerator and denominator by (pK 2 ), and use 
momentum conservation to eliminate p, for example 

[lp]<p3> = [l[p|3> = [1|2|3> = [12]<23> 
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this results in 

A(K+ K- v+ v+ v+) [12]3 1 (K2M + HM(Ki\t\K 2 ) 

which coincides indeed with the factorization limit of the exact expression calculated (gener- 
ally there are additional terms in the exact expression which are non-singular in the limit). 
The second factorization limit is for the channel Q = p% + K\ , in that limit 

A 5 (K 1 ,K 2 ,p+,p+,pt)^A 3 (pt,K 1 ,-Q)—A i (Q,K 2 , P +,p+) 

where Q is the momentum of a D-dimensional intermediate fermion. However physical 
amplitudes such as As(Ki, K 2 ,pf ,p 2 ,p%) are to be assembled from their components such as 
A&{Ki, K 2 , pt j p£ ) P3) which do not separately obey factorization constraints. We therefore 
did not check these more complicated factorization limits involving massive intermediate 
momentum. 
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